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^\J , Abstract. This work analyzes the distribution and size of interparticle gaps 

arising in an ensemble of hexagonal unit structures in the xy plane when pack- 
ing disks with a Gaussian distribution of radii with mean (r) and standard 
i-S^ ' deviation Ar. During the course of this investigation an equivalency is estab- 

lished between gaps arising in hexagonal unit structure packs and nine-ball 
a^^ ■ billiard rack patterns. An analytic expression is derived for the probability 

' distribution and location of interparticle gaps of magnitude F. Due to the 

number of variables and large number of possible arrangements, a Monte Carlo 

jrt ■ simulation has been conducted to complement and probe the analytic form for 

■4-J 1 three very different systems: i) billiard balls with Billiard Congress of America 

, ' (BCA) specifications, ii) US pennies with specifications of the US Mint, and 

iii) a hypothetical system with r = 1.0m and Ar = IzlO"^" m corresponding 

e^^ , to the scale of one atomic radius. In each case, probability density distribu- 

■ tions of gap sizes have been calculated for those Ar above, and also for 2Ar 

I ' and 0.5Ar, respectively. A general result is presented for the probability of a 

^ ' nonzero normalized (3-) gap size arising, P(r > oAr) = 1 — 0.124o, where a 

Sh , is a constant < 5.0. This curious result reflects the phenomenon of geometric 

w . frustration; the inability of the system to simultaneously satisfy all geometric 

constraints required by a perfect-rack sans interparticle gaps. 



00 I 1. Introduction 

cn ■ 

—1. . 1.1. Motivation. This work analyzes the distribution and size of interparticle gaps 

arising, locally, in an ensemble of non-interacting, hexagonal unit structures in the 
xy plane when packing disks with a Gaussian distribution of radii. During the 
r — ' course of this investigation an equivalency is established between gaps arising in 

(^ . hexagonal unit structure packs and nine-ball billiard rack patterns, see Figure 1. In 

fact, the motivation behind this work was the curious difficulty in preparing a gap- 
free rack of nine-ball billiards. This study maintains that upon packing disks with 
a Gaussian variation in radii, a perfect-pack, hexagonal and planar, is an illusion 
and is not attainable in practice due to geometric frustration. A perfect rack being 
one in which each disk contacts each of its neighbors with no gaps present. This 
paper explores this curious observable phenomenon for the simplified case of a nine- 
ball rack of disks with radii that vary according to a Gaussian distribution with 
standard deviation, a = Ar. Variation in disk radii (Ar) arising naturally due to 
usage and/or limitations in mill tolerance. 

Due to the initial motivation for this study, gaps arising due to local variables (i.e. 
random variation in disk radii) have been the primary focus of this study. Thus, 



Date: November 5, 2007. 



2 DANIEL P. SNOWMAN 

the term, ensemble, refers to a large number of non-interacting unit structures as 
shown in Figure 4. 

Intriguing parallels exist between this investigation and other physical systems 
with frustration present. Frustration in this context refers to the inability of the 
system to simultaneously satisfy all constraints present. Typically, in physical sys- 
tems, these constraints are related directly to the energy and entropy of the system. 
The present study is subject to geometric frustration that arises in that three circles 
(see below for development), and their respective equations, must be satisfied by a 
single {x,y}. Ultimately, geometric frustration arises if the system with circles of 
three different radii and centers do not share a common intersection point. The cen- 
ter and radii of each of these circles having an element of randomness present. Thus, 
this study analytically explores the nature of this geometric frustration that reveals 
itself as an interparticle gap in ordered, hexagonal packings of disks in the plane. 
In particular, this investigation, i) derives an analytic expression for the interpar- 
ticle gap of magnitude F with an associated probability distribution, ii) explores 
all physically realizable configurations of hexagonal disk packs, detailing resulting 
interparticle gap locations (see Figure 5), and, iii) performs a Monte Carlo simu- 
lation to probe the resulting probability density distribution of interparticle gaps 
arising in an ensemble of 10^ units. 

1.2. Disk Packings. Physicists, chemists, engineers, geologists, biologists and 
mathematicians have found many applications of hard disk packings as simipli- 
fied models for the study of granular materials [3] , random media [H] , glasses [SI] , 
liquids [12], amorphous solids [24], cell pattern formation, investigation of island 
formation in metallic films [4] , segregation problems [23] , turbulence [16] and plate 
tectonics [13]. The reader is directed to Aste and Weaire [2] for a quality account 
of many challenging theoretical problems and applications related to the field of 
hard disk packings. 

An earlier study by Cowan [5] employed a random packing strategy, following 
the Renyi parking scheme, as equal-sized disks were placed along the circumference 
of a central disk. This investigation probabilistically probed the number of gaps 
arising when the central disk was in contact with K neighbors. Cowan's study, 
much like the work presented in this paper, limited its scope to the immediate 
neighborhood of any one disk of the larger ensemble. Cowan later [7] considered 
the entire ensemble and established bounds on the density of the disk packing. 
The study presented in this paper differs as we consider ordered packing of disks 
with a Gaussian distribution of radii. The randomness lies in the disk-size, not 
the packing strategy. Aste [1] considered disordered packing of disks with random 
sizes and positions with the primary result being a power-law distribution for the 
size distribution of disks that densely cover the plane. Most intriguing in this 
study is the agreement found with an experimental investigation involving the vapor 
deposition of tin drops onto a hot substrate resulting in breath figures. 

Many of the more recent studies in the field of hard disk or sphere packing have 
dealt with the phenomenon of jamming [52], [IH]. In most cases, these packings 
consider packing fractions [S], density fluctuations [5] and possible glass transi- 
tions [9] for various classification of jammings [22], for a large number of particles. 
Our study does not directly investigate jamming as we consider a large number 
of noninteracting unit structures of ordered, hexagonal disk packs with a random 
distribution of disk radii. The results of the present study, however, may find direct 
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Figure 1 . A traditional nine-ball rack pattern. 

application to these works as the distribution of interparticle gap sizes in these unit 
structures directly affect the response of granular packings when subjected to a 
disturbing force. , , 

1.3. Circle Packings. A conjecture by Thurston [50" proposed an intimate con- 
nection between circle-packings and analytic function theory; that is, the rigidity 
associated with analytic functions also exists with circle-packings. Rodin and Sul- 
livan proved the corresponding convergence [17j . Many other profound results have 
been forthcoming since and the interested reader is directed to Stephenson for a 
comprehensive [I9] review of the field, or, see [H] for a short introduction to the 
topic. 

It is important to note that a distinction exists between a circle packing and 
sphere packing in two dimensions. Circle packings have an underlying pattern of 
tangencies, essentially consisting of an abstract graph, with a collection of vertices 
and connecting edges, that ultimately determines the circle packing pattern that 
arises as each circle adjusts its radius until the underlying pattern of tangencies is 
satisfied. 

The work presented in the present paper, technically, is not a circle packing 
since an underlying pattern of tangencies was not used to generate the packing. The 
reader should note, however, that after the ordered packing of randomly sized circles 
(or disks) is complete, a pattern of tangencies could be prescribed ex post facto. 
Thus, this study could be considered a pseudo-circle-packing investigation with 
randomness introduced in the underlying pattern of tangencies. This randomness 
would affect both features of the underlying associated graph: the length of each 
edge and the location of each node. 

2. Packing Strategy 

2.1. Nine-ball racks and packing order. Figure 1 shows a classic nine-ball 
diamond pattern rack of billiard balls. Only two balls (the one-ball and nine-ball, 
respectively) have predetermined positions in a proper nine-ball billiard rack; the 
one-ball at the head of the rack and the nine-ball at the center. The remaining 
seven balls, with enumerations 2-8, are placed in any of the remaining positions. 
It is interesting to note that a proper nine-ball rack can be prepared in 7! unique 
ways for one set of radii. Each unique configuration of disks potentially producing 
a unique gap size. 
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Figure 2. Most common interparticle gap locations in a nine-ball rack 

The billiard game of nine-ball demands the lowest numbered ball on the table 
always be struck first. Therefore, the one-ball would be located at the top of rack 
and typically is the first to be positioned. This assumes, of course, that the initial 
impact from the cue ball would be coming from the top of the page. After the one- 
ball is placed, a row of two-balls (with any enumeration 2-8) are placed snuggly 
behind the one-ball. A row of three balls with the nine-ball at the center is placed 
leaving no gaps. The placement of the fourth row of balls inevitably results in 
an interparticle gap arising in any one of the three positions shown in Figure 2. 
Finally, the remaining ball is placed. Our study involving disk packs in this same 
diamond shape will involve a slightly different packing strategy. 

2.2. Disk-packs and packing order. This study seeks to characterize the result- 
ing gap defects arising, thus, without any loss of generality we eliminate the disks 
at the top and bottom of the nine-ball pack shown in Figure 1, since neither con- 
tributes to gaps arising in our particular study. We are now left with a hexagonal 
structure as shown in Figure 3a and 3b. In our study, for computational ease, the 
nine-ball is placed first (at the origin) and from hereon is considered the center ball 
at the zeroth site as we transition away from the language of nine-ball billiards and 
into the arena of disk packing. 

2.3. Packing Strategy. After the center disk (or nine ball) is placed at the origin, 
the surrounding disks are packed in a counter-clockwise manner beginning at site 
si, as shown in Figure 3. The disk at site si lies directly on the x-axis a distance 
ri -I- To from the origin. The disk at site S2 is placed so that it is in contact with 
the disks at sites sq and si. Each pair of disks in contact share a common tangent. 
The disk at site S3 is placed so that it is in contact with the disks at sites so and S2. 
This process is repeated for the disks to be placed at sites S4 and S5 . It is assumed 
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Figure 3. Nine-ball rack pattern and equivalent hexagonal unit 
structure used in this study. 



that each disk to be placed does not disturb any disk previously placed. If the seven 
disks have identical radii, as shown in Figure 3b, a perfect rack will result with 12 
points of contact after the sixth and final disk has been placed. However, with disk 
of varying radii the final disk to be placed poses potential problems since the radii 
are not identical and it may not be possible to satisfy all geometric constraints 
present, given by Equations 5-7. 

3. Determination of Disk Location 

The location of each disk's center is determined by its radius, disk 0, and the 
disk placed immediately prior. The Law of Cosines allows for the determination of 
the angle that the i*'* disk's center will lie upon at a distance tq + Vi for each disk i. 
Figure 4b is drawn using the fact that each pair of disks in contact will have their 
two centers separated by a straight line of length r^ + tj . 

(1) % = arccos[(^'+'^^-^'7('^-+:°^'-(7+^°)'] 

^ ' y I -2(r, + ro)(r, +ro) ^ 

where i and j are immediate neighbors along the perimeter. Using this expression 
and Figure 3, we can write 



(2) ^15 = e 



12 + 1^23 -|- P34 + C45 
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Figure 4. Qualitative depiction of circle-of-contact equations rep- 
resenting the curve along which disks at sites 0, 1 and 5, respec- 
tively; will be in contact with a disk to be placed at site 6. 

Thus, in polar coordinates the center of disk 5 is located a distance r^ + rg from 
the origin at an angle ^15. 



(3) 

(4) 



X5 = (?'5 + ?'o) cos 6*15 

2/5 = {n +ro) sin 6*15 



The relationship derived above for {2:5, y^} is exact for one set of radii, {r^}, and 
single unit structure. However, a large ensemble of these unit structures, each with 
a unique set of radii, r^, requires the introduction of a probability distribution; see 
Section 5. 

After disks have been placed at sites 1 to 5 the final task remains: to place a 
disk at site six and determine analytically any resultant interparticle gaps arising. 
The placement of the disk at sites 2 to 5 is determined by the center disk and the 
disk placed immediately prior. The sixth and final disk is unique in that a perfect 
rack would have it contact three previously placed disk at sites sq, S5 and si. 

Contact between neighboring disks implies a separation between disk centers 
equal to the sum of their respective radii. Three circle-of-contact equations can 
be written with knowledge of their centers and corresponding radii. Equation 5 
represents a circle-of-contact between disks at sites and 6 with the center located 
at {xq^tjq} = {0,0} with radius rg + rg. Second, equation 6 represents a circle- 
of-contact between disks at sites 5 and 6 with center located at {x^jij^} = {(rg -I- 
r5)cos0i5, (tq -\- rs) sin 6*15} with radius r^ + rg. Finally, equation 7 represents a 
circle-of-contact between disks at sites 1 and 6 with center located at {xi,yi} = 
{fi -I- ro, 0} with radius ri -I- rg. 

The circle-of-contact equations are used to analytically determine all possible 
placements of one disk relative to its neighbor that results in contact between the 
two. Visually, circles-of-contact also have a strong appeal, as shown in Figure 4. 
Figure 5 qualitatively depicts all types of intersections occurring between the three 
circle-of-contact equations. 

The location of the sixth disk {xQ,yQ} must satisfy each of the three circle- 
of-contact equations. Equations 5-7, if it is to be in contact with all three of its 
neighbors at sites 0, 5, and 1. A perfect rack will result only if these three circle- 
of-contact equations share a common intersection point. If only two of the three 
relations intersect an interparticle gap will be present. 
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(5) x^+y^ = (ro + Vf^f 

(6) [x - x^f + {y - y^f = {rr-, + rf,f 

(7) {x-xif+y^ = {ri+rQf 

This study is concerned with finding those {x, y} that satisfy at least two of the 
three equations above. A triple intersection point would correspond to an {x, y} 
that satisfies all three of the equations, corresponding to a perfect rack with no gap 
defect. 

4. Possible Configurations and Gap Locations 

Figure 5 depicts twenty-five nontrivial, at least two of the three circle-of-contact 
equations share a solution, intersections in the fourth quadrant. The fourth quad- 
rant is the focus of this study as we attempt to place the sixth disk, see Figures 
3 and 4, subject to the functional constraints of three-separate circle-of-contact 
equations given by Equations 5-7. Those intersection points, shown in Figure 5, 
represent the location of the center of disk six resulting in contact with two of its 
three neighboring disks. 

This investigation further reduces the number of solutions considered as we elim- 
inate those intersections not physically possible. That is, if the location of any one 
of the intersection points is located inside any one of the three circle-of-contact 
curves, that solution is eliminated since it is not physically possible. That is, the 
closest approach of any two disk centers is a sum of their respective radii. Inspection 
of Figure 5 reveals that intersection points e, f, g, i, j, q, and y are not physically 
realizable since these intersections would demand that two neighboring disks are 
separated by a distance less than the sum of their respective radii. 

Of the 18 remaining solutions: i) three triple intersection points exist (p, v, and 
x) corresponding to zero interparticle gap, ii) five points of intersection (a, k, m, r 
and t) leading to a gap between sites 5 and 6, Fsg, iii) five points (b, 1, n, s, and 
u) of intersection leading to a gap between sites 1 and 6, Fig, and iv) five points of 
intersection (c, d, h, o and w) leading to a gap between sites and 6, Fog. Of these 
eighteen physically realizable configurations, only 13 are unique, with symmetry 
arguments allowing us to qualitatively equate the following; a=b, k=l, m=n, r=s 
and t=u. 

5. Analytic Gap Size 

After determining the disk configurations physically possible (see Section 4) for 
a given set of radii, the size of any interparticle gap present between disk six and 
any one of its neighbors is determined using Equations 8-10. 

(8) Foe = [{xa - x^f + (yo - yef] ^ - (ro + rg) 

(9) Fi6 = [(xi - x^f + (yi - y^f] 5 - (n + re) 

(10) Fse = [(xs - x^f + (yg - y^f] ^ - (rg + rg) 

Regardless of the configuration resulting, dependent upon the particular set of 
Ti, a minimum of two of the three Equations 8-10 will be zero, corresponding to 
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Figure 5 . Qualitative depiction of all possible intersections of the 
three circle-of-contact equations. Each intersection represents the 
potential, future location of disk 6 

those neighboring disks in contact with disk six. The nonzero F represents the 
interparticle gap size occurring with probability P{T), shown in Equation 11. 



(11) 

with 
(12) 
and 

(13) 



P(F) = l[P{r 



i=0 



a — Ar 



Pin) - ^i?xp[^^-5^ 



27rcr 
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6. Monte Carlo Analysis 



The probability distribution for the gap size depends upon each of the seven radii, 
thus, due to the number of variables and the large number of possible arrangements, 
a numerical simulation has been conducted to complement and probe the analytic 
form of the gap probability distribution. 

Each simulation involves generating seven radii from a Gaussian distribution 
and packing them counter-clockwise starting at site 1 (as described in Section 2). 
For each of 10^ packings, a gap size was calculated and a resultant probability 
density distribution produced. This study has considered three very different disk- 
packings with disks of radii that follow a Gaussian distribution with mean (r) 
and standard deviation (Ar). The systems considered were: i) billiard balls with 
Billiard Congress of America (BCA) specifications, ii) US pennies with specification 
of the US Mint, and iii) a hypothetical system with r = 1.0m and Ar = lxlO~^'^m 
corresponding to the scale of one atomic radius. In each case the probability density 
distribution of gap sizes has been calculated for those Ar above, and also for 2Ar 
and 0.5Ar. 

7. Results 

In each case detailed below, it was found that the distribution of gaps arising 
between disk 6 and neighbors 0, 1, and 5, respectively (and computationally calcu- 
lated with Foe, Fig, and Fsg) are statistically indistinguishable from one another. 
Thus, the analysis below pertains to gap size distributions regardless of location. 

For each of the three systems considered, the resulting probability density distri- 
butions were found to be Gaussian in nature. Normalization of these distributions 
for each of these unique systems demands the integral shown in Equation 14 equal 
1 over the range of positive F. In each case, the best-fit Gaussians were compared 
to the distributions resulting from the Monte Carlo simulation via a qqplot. Corre- 
lation coefficients were calculated for each case and showed exceptional agreement 
with the Gaussian fits, see Section 7.4 and Figure 14. 



(14) / pExp[-^^]dr^l 

Jo 



2a^' 



Evaluating the integral of Equation 14 leads to an expression for the normaliza- 
tion constant 

(15) (3- 



/2TTa 

Thus, each probability density distribution presented below has a Gaussian of 
the form shown in Equation 16. 

(16) PDF{r, a) = -^Exp[-^] 

\/2tt(t 2(7^ 

7.1. Billiard Balls. Figure 6 represents the probability density as a function of 
gap size for a disk packing with disk radii and mill tolerances identical to those 
specified by the Billiard Congress of America (BCA); that is, billiard balls used in 
competition must have radu 28.57500 ± 0.06350TOm. 
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Figure 6. Probability density distribution for interparticle gap 
sizes arising in a packing of billiard balls meeting Billiard Con- 
gress of America (BCA) specifications, r — 28.57500?tt-77i and 
Ar = 0.06350mTO, with best-fit Gaussian. 



The probability density data shown in Figure 6 is described exceptionally well 
with a Gaussian best-fit with the functional form of Equation 16. The proba- 
bility density data plotted and the best-fit curve both have standard deviation 
a — 1.63a;10~^ mm. Figure 7 compares gap size probability density distribu- 
tions arising using the BCA specification Ar = 0.06350TO?Ti to those resulting 
when: i) Ar is doubled to 0.12700TOm, and ii) Ar is halved to 0.03175?tito. Dou- 
bling Ar results in a broader distribution characterized with a standard deviation 
a = 3.25xlO~^?TiTO, consistent with both simulation data and best-fit Gaussian. 
Halving Ar results in a sharper gap size distribution characterized with a standard 
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Figure 7. A comparison of probability density distributions for 
interparticle gap sizes arising in packing of billiard balls with r — 
28.5750mm and i) doubled BCA spec Ar = 0.12700?ti?tt, denoted 
by (a), ii) BCA spec. Ar — O.OGSSOmrn denoted by (•), and iii) 
halved BCA spec, to Ar = 0.03157toto denoted by (■); with 
best-fit Gaussians. 



deviation a = 8.18x10 



consistent with both Monte Carlo simulation data 



and best-fit Gaussian. The reader should note that a doubling (halving) of Ar 
results in a proportional change in the standard deviation, a, characterizing the 
probability density distribution. 

The average gap size calculated analytically using the best-fit Gaussian and 



Equation 17, reveal Fq 
simulation data. 



1.29x10 mm, in agreement with the Monte Carlo 



12 DANIEL P. SNOWMAN 



(17) r^ve = / r^=^Exp[--^]dr 

Jo VSttct 2ct^ 

7.2. United States Pennies. Figure 8 represents the probability density as a 
function of gap size for a disk packing with disk radii and mill tolerances identical 
to those specified by the United States Mint for the manufacture and production 
of US Pennies; that is, US Pennies must have radii 9.5250mm ± O.OOSOmm. This 
may be the most interesting of the three cases considered in this paper since an 
actual, experimental disk-packing can be easily performed using seven pennies and 
a tabletop. Reader beware, however, the average interparticle gaps arising is shown 
below to be only 2.60xlQ~^mm. 

The probability density distribution shown in Figure 8 hereto is represented 
very well with a Gaussian best-fit with the functional form of Equation 16. The 
probability density data plotted and the best-fit curve both have standard deviation 
(T — 3.26xl0^^mm. Figure 9 compares gap size probability density distribution aris- 
ing using the US Mint specification Ar = 0.0050m?Ti to those resulting when: i) Ar 
is doubled to O.OlOOmm, and ii) Ar is halved to 0.0025TOm. Doubling Ar results in 
a broader distribution characterized with a standard deviation a = 6A5xlO~^mm, 
consistent with both simulation data and best-fit Gaussian. Halving Ar results in 
a sharper distribution characterized with a standard deviation a — 1.62xl0^^mm, 
consistent with both simulation data and best-fit Gaussian. Once again, notice that 
doubling (halving) Ar results in a corresponding doubling (halving) of the standard 
deviation, cr, the characteristic parameter defining the sharpness of the probability 
density distribution. 

The average gap size calculated analytically using the best-fit Gaussian and 
Equation 17, reveal Tave — 2.60xl0^^mm, in agreement with the Monte Carlo 
simulation data. 

7.3. Hula-Hoops and Angstroms. Figure 10 represents the probability density 
as a function of gap size for a disk packing with disk radii 1.0m ± 10~"'^''m. This 
particular disk-packing was included in this study as it would represent an extreme 
physical minimum for Ar; that is, the characteristic scale for atomic radii. This 
disk packing could be analogous to packing hula-hoops of radius l.Om with a very 
unrealistically tight mill tolerance of 1 Angstrom. 

As with the previous two studies discussed in Sections 7.1 and 7.2 above, the 
probability density distribution shown in Figure 10 is well represented with a Gauss- 
ian best-fit with the functional form of Equation 16. The probability density data 
plotted and the best-fit curve both have standard deviation a = 6.49a;10~^'' m. Fig- 
ure 1 1 compares gap size probability density distribution arising using the angstrom 
specification Ar = 1.0x10""'^'' m to those resulting when: i) Ar is doubled to 
2.0x10"^'' m, and ii) Ar is halved to 5.0a;10~"'^^TO. Doubling Ar results in a broader 
distribution characterized with a standard deviation a = 1.30a:10~^m-, consistent 
with both simulation data and best-fit Gaussian. Halving Ar results in a sharper 
distribution characterized with a standard deviation a = 3.23a;10~"'^''TO, consistent 
with both simulation data and best-fit Gaussian. Consistent with the billiard and 
penny studies above, notice that doubling (halving) Ar results in a correspond- 
ing doubling (halving) of the standard deviation, a, the characteristic parameter 
defining the sharpness of the probability density distribution. 
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Figure 8. Probability density distribution for interparticle gap 
sizes arising in a packing of pennies with US Mint specifications 
r — 9.5250TO?n and Ar ~ O.OOSOttito, with best-fit Gaussian. 



The average gap size calculated analytically using the best-fit Gaussian and 



Equation 17, reveal Ta 
simulation data. 



5.17x10 



-10^ 



in agreement with the Monte Carlo 



7.4. Universal Probability Density Distributions. To allow comparison of 
the three very different ranges and corresponding disk packings detailed above, a 
normalized probability density distribution was generated as a function of normal- 
ized gap size, ct = ^- The best-fit Gaussian, shown on Figure 12, also uses a 
normalized standard deviation a' = a j I\r and is shown as the smooth curve. The 
data for the three normalized probability density distributions exhibit exceptional 
agreement and all three of these probability density distributions are described by 
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Figure 9. A comparison of probability density distributions for 
interparticle gap sizes arising in packing of United States pennies 
with r — 9.525mTO and i) doubled US Mint spec Ar — O.OIOOtotti 
denoted by (a), ii) US Mint spec. Ar = O.OOSOttito denoted by 
(•), and iii) halved US Mint spec, to Ar = 0.0025toto denoted by 
(■); with best-fit Gaussians. 



the same best-fit Gaussian with a — 6.43. This figure reveals the universal nature 
of the geometric frustration that arises in these disk-packings regardless of scale. 
Figure 14 shows the qqplot for each of these normalized distributions comparing 
the distribution arising from the Monte Garlo simulation to the best-fit Gaussian. 
The result for this case yields a correlation coefficient of R = 0.999991 between the 
two distributions. Each of the previous distributions yields a correlation coefficient 
at least as good when analyzing each distribution using similar qqplot analysis. 
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Figure 10. Probability density distribution for interparticle gap 
sizes arising in a packing of disks with r = l.Om and Ar = 
1.0xlO^^"m, with best-fit Gaussian. 



The average gap size calculated analytically using the best-fit Gaussian and 
Equation 17, reveal Tave ~ 5.09q:, in agreement with the Monte Carlo simulation 
data. 

7.5. Probability of a gap size F > aAr. The probability of an interparticle 
gap arising, F > aAr, can be determined using standard probability theory, see 
Equation 18. 



(18) 



P{T > aAr) = 1- f -^Exp[^^]dT 
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Figure 11. A comparison of probability density distributions for 
interparticle gap sizes arising in packings of disks with r = 1.0 ni 
and i) Ar = 2.0x10~i°to denoted by (a), ii) Ar = LOxlO^^^m 
denoted by (•), and iii) Ar = 5.0x10~^^to denoted by (■); with 
best-fit Gaussians. 

This integral has been evaluated for a range of a, and the results are displayed 
in Figure 13. The probabilities shown represent those for the universal probability 
distribution with a = 6.43. The reader should note for a < 5.0, the probability of 
a gap arising F > aAr is approximately 1 — 0.1 24q; as determined by Figure 13. 
For increasing a (i.e. > 5.0) the probability becomes nonlinear as shown in Figure 
13. 



(19) P(r > aAr) = 1 



2 r^ 

Exp\ -]dT ' 



1- 0.124a /or a < 5.0 
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Figure 12. A comparison of normalized probability density dis- 
tributions as a function of normalized gap sizes for i) billiard balls 
with BCA specs, ii) US Pennies with the specs of the US Mint, and 
iii) an extreme disk packing with r = 1.0m and Ar — 1x10" 
with best-fit Gaussian. 



-10 



8. Conclusions 

8.1. General Comments. This work analyzes the distribution and size of gap de- 
fects arising, locally, in an ensemble of non-interacting, hexagonal unit structures in 
the xy plane when packing disks with a Gaussian distribution of radii. An analytic 
expression for the probability distribution and location of the gap defect of magni- 
tude r has been derived. Due to the complexity of the expression a Monte Carlo 
simulation was employed to probe the resulting frequency and size distribution of 
gaps arising in an ensemble of 10^ unit structures. A best-fit of the Monte-Carlo 
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Figure 13. Probability of an interparticle gap arising with mag- 
nitude r > aAr, for a range of a with fixed Ar. 



produced frequency distributions were performed and, in general, are represented 
very well with Gaussian best-fits as shown by a qqplot analysis for each case. 



8.2. Billiards and Granular Applications. Any billiards aficionado realizes the 
importance of the break in any game of pool, especially nine-ball. A novice will 
seek a re-rack to eliminate any gap that has suddenly appeared in a rack about to 
be scattered. The experienced player realizes the opportunity the gap represents. 
The presence of interparticle gaps drastically affect resultant scattering. General 
trends can be detected for the initial direction of most of the balls in a nine-ball 
rack immediately after the break, primarily from experience. However, an analytic 
expression accompanied with a probability distribution is the desired result. 
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Figure 14. A qqplot comparing the actual distribution and the 
best-fit Gaussian for the case of normahzed probability density 
distribution was generated as a function of normalized gap size. 
Comparing the two distributions yields a correlation coefficient of 
R = 0.999991. 



The results of the present paper will be used in future studies as the propagation 
of shock waves through nine-ball billiard racks and hexagonal granular packs with 
interparticle gaps is considered. More generally, upon the application of a force the 
scattering response of tightly packed granular materials and the affect of inevitable 
interparticle gaps upon this scattering will be considered. 

Further, the results presented in this study were motivated by curiosities arising 
in billiard ball racks, however, the replacement of spheres with disks is a simplifi- 
cation to be removed in an upcoming study. That is, an upcoming paper addresses 
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gap defects arising in a packing of spheres, with a Gaussian distribution of radii, 
on the plane. 

8.3. Disk Packing Application. The results of this work may also have applica- 
tion to Cowan type disk packings. Whereas Cowan considered a random packing 
strategy for identical disks, the present study considered ordered packing of disks 
with random elements affecting the disk-size, not the packing strategy. 

Second, the probability distributions presented in this paper for interparticle 
gap size may find use in studies considering polydisperse disk packings covering the 
plane. 

In addition, although this work has not been an investigation of jamming in disk 
packs, the results may find application in this field since the interparticle gap size 
in these unit structures may directly affect the response of these packings upon the 
application of a disturbing force. 

8.4. Circle Packing Application. The results of this work may have application 
to field of circle packing if one were to consider the possibility of an underlying 
pattern of tangencies with random elements affecting the length of each edge and 
location of each node. 
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